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ABSTRACT 

Supermassive primordial stars are expected to form in a small fraction of massive protogalaxies 
in the early universe, and are generally conceived of as the progenitors of the seeds of supermassive 
black holes (BHs) at high redshift. Supermassive stars with masses of ~ 55,000 M©, however, have 
been found to explode and completely disrupt in a supernova (SN) with an energy of up to ^ 10^^ 
erg, instead of cohapsing to a BH. Such events, roughly 10,000 times more energetic than typical 
SNe today, would be among the biggest explosions in the history of the universe. We carry out a 
simulation of such a supermassive star SN in two stages. Using the RAGE radiation hydrodynamics 
code we first evolve the explosion from the earliest stages, through the breakout of the shock from 
the surface of the star until the blast wave has propagated out to several parsecs from the explosion 
site, which lies deep within an atomic cooling dark matter (DM) halo at z ~ 15. Then, using the 
GADGET cosmological hydrodynamics code we evolve the explosion out to several kiloparsecs from 
the explosion site, far into the low-density intergalactic medium. The host DM halo, with a total 
mass of 4 X 10 Mq, much more massive than typical primordial star- forming halos, is completely 
evacuated of high density gas after < 10 Myr, although dense metal-enriched gas recoUapses into the 
halo, where it will likely form second-generation stars after > 70 Myr. The ^ 20,000 Mq in metals 
that are released in the explosion are widely distributed, and enrich the dense recoUapsing gas to an 
average metallicity of ~ 0.05 Zq. Such a high level of enrichment suggests that the chemical signature 
of these supermassive star explosions may have been missed in previous surveys of metal-poor stars. 
Subject headings: Cosmology: theory — early universe — supernovae: general 



1. INTRODUCTION 

Recently, there has been renewed interest in the long- 
standing theoretical possibility that supermassive stars 
(SMSs), with masses of 10^-10^ Mq (e.g., Iben 1963; 
Fowler & Hoyle 1964; Appenzeller & Fricke 1972; Shapiro 
& Teukolsky 1979; Bond et al. 1984; Fuller et al. 1986), 
inhabited the early universe (see e.g., Volonteri 2012). 

One of the main motivations for their study comes from 
observations of quasars at z ~ 6-7 which are inferred to 
be powered by black holes (BHs) with masses exceed- 
ing 10^ Mq (e.g., Willott et al. 2003; Fan et al. 2006; 
Mortlock et al. 2011). Given the short time (< 800 Myr) 
available for such massive BHs to grow via accretion from 
their initial 'seed' masses, as derived from the most recent 
cosmological parameters inferred by e.g., the Wilkinson 
Microwave Anisotropy Probe (Komatsu et al. 2011)[3 
and the suppression of BH growth due to the strong 
radiative feedback from both stars (e.g., Whalen et al. 
2004; Wise & Abel 2007; O'Shea & Norman 2008) and 
the BHs themselves (e.g., Pelupessy et al. 2007; Alvarez 
et al. 2009; Milosavljevic et al. 2009; Jeon et al. 2012; 
Park & Ricotti 2012), it now appears more likely than 



ever that the seeds of the most massive early BHs must 
have been quite massive themselves (e.g., > 10^ Mq] see 
Johnson et al. 2012a; also e.g., Shapiro 2005; Volonteri 
& Rees 2006; Natarajan & Volonteri 2012). Whereas 
the majority of the first, Population (Pop) HI stars may 
have had masses of -- 20-500 Mq (e.g., Abel et al. 2002; 
Bromm & Larson 2004; Yoshida ct al. 2008; Greif et 
al. 2011), the best candidates for the seeds of SMBHs 
are thus much more massive (and rare) supermassive pri- 
mordial stars. 

An additional, and independent, reason to consider 
SMSs in the early universe is that the conditions re- 
quired for their formation are now thought to be real- 
ized much more often than was previously assumed. The 
most widely discussed avenue for the formation of SMSs 
is via the direct gravitational collapse of hot (~ 10'' K) 
primordial gas in so-called atomic cooling dark matter 
(DM) halos at z > 10 (e.g., Bromm & Loeb 2003; Begel- 
man et al. 2006; Lodato & Natarajan 2006; Spaans & 
Silk 2006; Regan & Haehneh 2009; Choi et al. 2013; 
Latif et al. 2013a, b)Q In this scenario, the gas in the 
protogalaxy remains at the virial temperature of ^ 10* 



^ Adopting the cosmological parameters reported recently by the 
Planck Collaboration (2013) yields a similar time available for seed 
growth. 



^ We note that other formation mechanisms stemming from 
shocks (Inayoshi et al. 2012) and magnetic fields (Sethi et al. 2010) 
have also been suggested. 
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K because H2 molecules have been photodissociated by 
the LW background, leading to the rapid formation of 
SMSs via the accretion of gas at rates ^ 10^-10'^ times 
higher than in the formation of most Pop III stars from 
H2-cooled gas. The flux of radiation required to keep the 
gas H2-free depends on its spectrum, with lower fluxes 
required if it is produced by metal-enriched stars instead 
of Pop III stars (e.g., Shang et al. 2010). Recent work by 
independent groups has shown that Pop II star-forming 
galaxies in the early universe are able to produce suffi- 
cient H2-dissociating radiation to prevent the cooling of 
primordial gas in a substantial fraction of atomic cool- 
ing halos, thereby leading to the seeding of these halos 
with SMSs that can collapse into BHs (see Dijkstra et al. 
2008; Agarwal et al. 2012; Petri et al. 2012; Johnson et 
al. 2013). Indeed, Agarwal et al. (2012, 2013) find that 
a large fraction of the SMBHs in the centers of galaxies 
today may have been seeded by SMSs. Strengthening 
these conclusions are other recent results which suggest 
that lower LW fluxes may be required for the forma- 
tion of SMSs, due to a reduced role of H2 self-shielding 
(Wolcott- Green et al. 2011) and the presence of signifi- 
cant turbulence or magnetic fields (Van Borm & Spaans 
2013). 

Complementary studies have been undertaken to un- 
derstand the growth and evolution of SMSs, as well. 
Modeling the growth of accreting protostars with masses 
up to ~ 10^ -^'^o, Hosokawa et al. (2012) have shown 
that they emit little high energy radiation that could halt 
their continued accretion, and Inayoshi et al. (2013) have 
shown that pulsational instabilities are likewise unable to 
halt their growth. Johnson et al. (2012b) modeled the 
growth of SMSs to much higher masses and showed that, 
even if they are able to emit the copious ionizing radia- 
tion characteristic of main sequence Pop III stars, radia- 
tive feedback is not able to stop their growth up to at 
least ~ 10^ Mq. At the highest accretion rates expected 
for these objects (> 1 Mq yr~-^; Wise et al. 2008; Shang 
et al. 2010; Johnson et al. 2011), the masses of primor- 
dial SMSs are only limited by the < 4 Myr that they 
have to accrete gas before they collapse to BHs (Begel- 
man 2010). 

Whereas the majority of SMSs are expected to collapse 
to black holes with little or no associated explosion (e.g.. 
Fryer & Heger 2011), it is possible that some fraction 
instead explode as extremely energetic supernovae (SNe; 
e.g.. Fuller et al. 1986; Montero et al. 2012; Whalen et 
al. 2012; see also Fuller & Shi 1998; Linke et al. 2001). 
In particular, Heger et al. (2013) have found that SMSs 
with masses in a narrow range around ~ 55,000 A/© end 
their lives as extraordinarily luminous SNe|j With ener- 
gies of almost 10^^ erg, these explosions are among the 
most energetic in the history of the universe. 

Here we expand on the radiation hydrodynamics simu- 
lations presented by Whalen et al. (2012a) and simulate 
the long-term evolution of a SMS SN in its cosmologi- 
cal environment, in order to show how these gargantuan 
explosions impact both the formation of the first galax- 
ies and the chemical signature of the first stars. In the 
next section, we describe the multi-scale simulations that 

^ These SNc would appear much brighter than other types of 
Pop III SNe (e.g., Scannapieco et al. 2005; Hummel et al. 2012; 
Pan et al. 2012; Whalen et al. 2012b,c,d). 
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Fig. 1. — The initial conditions for the GADGET 3-D cosmolog- 
ical simulation. The blue curves show the radial velocity {left) and 
hydrogen number density (right) profiles resulting from the 1-D 
RAGE simulation of the early evolution of the SMS SN remnant, 
up to the time when the shock extends out to 6 pc from the explo- 
sion site. The orange triangles denote the GADGET SPH particles 
representing the SN ejecta that are fit to the RAGE output. The 
black points denote the the SPH particles within the cosmological 
halo hosting the explosion. 



we have carried out to model the evolution of the explo- 
sion from its breakout from the surface of the star to 
the propagation of the blast wave into the intergalactic 
medium (IGM). In Section 3, we present our results on 
the energetics and dynamics of the explosion, as well as 
on metal enrichment and second-generation star forma- 
tion. In Section 4, we conclude with a brief discussion of 
our results. 

2. SIMULATION SETUP 

Here we describe the two simulations that we have 
carried out. The first is a 1-D radiation hydrodynam- 
ics calculation using the Los Alamos National Labora- 
tory RAGE code (Gittings et al. 2008) which allows 
to track the propagation of the SN blast wave out to 
several parsecs from the explosion site, deep within the 
host atomic cooling halo. For the second, we map the 
results of the first into a 3-D cosmological simulation us- 
ing the GADGET hydrodynamics code (Springel et al. 
2001; Springel & Hernquist 2002). Our use of these two 
simulation codes for the phases of the SN in which they 
are most well-suited to accurately model the explosion, 
from small (AU) scales to large (kpc) scales, constitutes a 
significant improvement over previous cosmological sim- 
ulations of SN feedback. 

2.1. Stellar and Early Supernova Evolution 

For the SMS progenitor of the SN we adopt the 55,000 
Mq stellar model described in Whalen et al. (2012a), 
which was evolved until the onset of explosion, using 
the Kepler code (Weaver et al. 1978; Woosley et al. 
2002). The explosion was then followed using Kepler 
and was confirmed using the CASTRO code (Almgren 
et al. 2010). The explosion completely disrupts the star 
and yields an explosion energy of 7.7 x 10^* erg (Heger 
et al. 2013). The RAGE code is then used to simulate 
the SN from the breakout of the shock from the surface 
of the star until the blast wave has propagated through 
a circumstellar medium with a density > 10^ cm~'^ (and 
a density profile ex r~^) out to several parsecs from the 
explosion site. Up to this point, the stellar evolution cal- 
culation and the simulation of the early phases of the SN 
are the same as described in Whalen et al. (2012a), to 
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Fig. 2. — Properties of the gas in the vicinity of the SN at 100 kyr (top panels), 1 Myr [middle panels) and 10 Myr {bottom panels) 
after the explosion of the central SMS. From left to right, as functions of the distance from the explosion site, are the number density of 
hydrogen nuclei, outward radial velocity, gas temperature and free electron fraction. Contours denote the distribution of the gas, with the 
mass fraction varying by an order of magnitude across contour lines. The ~ 1000 km s~^ shock completely disperses the dense (n > 10 
cm"'^) gas in the center of the host atomic cooling halo, but is only slightly decelerated even after propagating well beyond the ~ 1 kpc 
virial radius of the halo. 



which we refer the reader for more detailed discussion of 
the calclulations (see also Frey et al. 2013). 

The blue curves in Figure 1 show the velocity (left 
panel) and density (right panel) profiles of the SN ejecta 
at the end of the RAGE simulation, at which point the 
shock has propagated out to 6 pc from the explosion site. 
The 55,000 Mq of ejecta, 23,000 Mq of which is heavy 
elements produced during the evolution and explosion of 
the progenitor, are traveling at almost 10,000 km s~^ 
outward from the center of the host atomic cooling halo. 
It is from this point that we map these velocity and den- 
sity profiles into a self-consistently evolved atomic cool- 
ing halo in a much larger cosmological volume. 

2.2. Cosmological Blast Wave 

To simulate the subsequent evolution of the SN blast 
wave in the appropriate cosmological environment, we 
map the velocity and density profiles obtained from the 
smaller-scale RAGE simulation into the center of a 4 x 
10^ Mq atomic cooling DM halo, the type of which is 
expected to host the formation of SMSs in the early uni- 
verse. The halo is identified in a 1 Mpc'^ (comoving) 



cosmological volume which has been evolved from z — 
100 down to z ~ 15 under the influence of a uniform, 
elevated H2-dissociating (Lyman- Werner; LW) radiation 
field, which prevents the gas from cooling and is assumed 
to lead to formation of a single SMS. Further details of 
the cosmlogical simulation up to this point are described 
in Johnson et al. (2011), who considered the impact of 
the alternative end state of such a SMS, a rapidly accret- 
ing BH. 

To map the output of the (Eulerian) RAGE simulation 
into the (Lagrangian) smoothed particle hydrodynamics 
(SPH) GADGET simulation, we assigned the central 460 
SPH particles, constituting the 55,000 AIq of gas within 
~ 6 pc of the densest particle in the halo, outward (ra- 
dial) velocities and hydrogen number densities so as to 
match those from the RAGE simulation. The fits that 
we obtain are shown in Fig. 1, with the inner 55,000 
Mq in SPH particles constituting the ejecta denoted by 
orange triangles and the unpertured particles residing in 
the outskirts of the halo denoted by black circles. While 
we fit the velocity profile very well, due to the mapping 
from an Eulerian to Lagrangian code the density profile is 
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Fig. 3. — The mass of material swept up by the shock {red dashed 
line) and the total kinetic energy of the swept-up material {green 
solid line), as a function of time since the SN. The majority of the 
kinetic energy is radiated away within the first 10* yr, and after 
50 Myr all laut ~ 1 percent has been lost. The majority of the 
material that is overtaken by the shock is swept-up after 1 Myr, 
with the total swept-up mass approaching 10^ Mq. 



somewhat noisiero Nonetheless, the basic features of the 
density profile containing the vast majority of the mass 
are represented, and the overall energy and momentum 
are also well-matched. From the initial conditions shown 
in Fig. 1 we restart the cosmological SPH simulation, 
the results of which we present in the next section. 

We note that, beyond mapping into it the blast profile 
of the SMS SN, we have chosen to leave the gas in the 
host halo otherwise unchanged. We have made this sim- 
plifying choice, in light of the large uncertainties in the 
radiative output of rapidly accreting SMSs, due to which 
it is unclear how the radiation emitted during the brief {r^ 
2 Myr) lifetime of the star will impact the medium within 
the host halo. Even if the star emits copious ionizing ra- 
diation, as main sequence Pop III stars are expected to 
do, it may be that the H ii region created by the star 
is confined to the innermost regions of the halo (John- 
son et al. 2012; see also Hosokawa et al. 2012 on the 
possibility of even less energetic radiation being emitted 
during the protostellar phase). It is possible, on the other 
hand, that ionizing radiation is able to escape out into 
the halo, if the accretion flow is highly anisotropic (e.g., 
due to the presence of an accretion disk; e.g., McKee & 
Tan 2008), or if accretion is intermittent (e.g., Clark et 
al. 2011; Smith et al. 2011; Vorobyov et al. 2013), in 
which case radiation could break out during periods of 
reduced accretion. 

Whereas the small-scale structure of the interstellar 
gas in the simulation we present here is subject to the 
limited resolution of the cosmological simulation into 
which we place the expanding blast wave, we note that 
at sub-resolution scales (< 1 pc) it is possible that a 
substantial amount of energy in the explosion is radi- 
ated away (e.g., Kitayama & Yoshida 2005; Whalen et 

* Note in Fig. 1 that the density spike at the shock front in 
the RAGE output contains less mass than is contained in a single 
GADGET SPH particle, and so there are no particles representing 
this particular parcel of gas in the cosmological simulation. This 
illustrates the fundamental difficulty in matching output from an 
Eulerian, adaptive mesh refinement code to an SPH code. 



al. 2008). We shall address how such additional radia- 
tive losses would affect the dynamics of the blast wave 
and the metal enrichment of the host halo and IGM in 
future work. 

3. RESULTS 

Here we present the results of our cosmological hydro- 
dynamics simulation, with particular attention paid to 
the dynamics of the expanding blast wave and to the 
enrichment of the IGM by the metal-rich SN ejecta. 

3.1. Dynamics and Energetics 

The injection of almost 10^^ erg at a single explosion 
site has a dramatic impact on the host halo. Figure 2 
shows the properties of the gas in the vicinity of the 
halo, as a function of the distance from the explosion 
site at its center. As the left-most panels show, the blast 
results in the complete evacuation of high-density (n > 
10 cm~^) gas within 10 Myr, with the material overtaken 
by the shock being carried out beyond the virial radius 
of the halo (at r ~ 10'^ pc) at up to ~ 10^^ km s^^. This 
material is also shock heated to temperatures up to ~ 10* 
K, resulting in its almost complete ionization, as shown in 
the right-panels. The gas, however, rapidly cools due to 
inverse Compton scattering of CMB photons, H and He 
atomic line emission, and bremsstrahlung, as also found 
(using the same GADGET code) in the less-energetic 
(10^^ erg) pair instability supernova (PSN) explosion in 
a 2.5 X 10^ DM halo simulated by Greif et al. (2007). 

As shown in Figure 3, most (~ 90 percent) of the 7.7 x 
10^^ erg initially in the blast is radiated away via these 
processes within 10^ yr. Nevertheless, the momentum of 
the blast is conserved and the shock continues to propa- 
gate into the IGM, sweeping up the majority of the mass 
after ~ 1 Myr. By 50 Myr, at which time ~ 99 percent of 
the energy has been radiated away, the shock has prop- 
agated out to ~ 5-10 kpc and has swept up almost 10^ 

Figure 4 shows the properties of the gas in the vicinity 
of the explosion site within a 400 pc (comoving) slice of 
the cosmological volume, at 1 Myr, 10 Myr and 50 Myr 
after the explosion of the SMS. Comparing the radial ve- 
locity field (second column from the left) to the cosmo- 
logical density field (far left column), it is clear that the 
blast wave propagates most rapidly into the low-density 
voids while its progress is halted in the direction of the 
high-density filaments, at the intersection of which lies 
the host halo. Figs. 2 and 4 also show the same general 
trend that, at late times, the most strongly shock- heated 
and highest- velocity material is located behind the shock 
front several kpc from the explosion site. The cooler ma- 
terial within the (~ 1 kpc) virial radius of the host halo 
is able to begin rccoUapsing after 50 Myr. As wc discuss 
next, this gas is likely to form second-generation stars 
that are enriched to fairly high metallicities. 

3.2. Metal Enrichment and Second- Generation Star 
Formation 

As is also expected for less energetic PSNe from Pop III 
stars with masses of ~ 200 Mq (e.g., Heger et al. 2003), 
a large fraction of the ejecta from our SMS SN consists 
of heavy elements. In particular, ~ 23,000 AIq of newly- 
synthesized metals are ejected in the explosion. These 
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Fig. 4. — Properties of the gas in the vicinity of the SN at 1 Myr [top panels), 10 Myr (middle panels) and 50 Myr (bottom panels) after 
the explosion. From left to right, are projections of the number density of hydrogen nuclei, outward radial velocity, gas temperature and 
free electron fraction. The lighter colors denote higher values for the quantities shown (see also Fig. 2). The explosion is asymmetric, with 
the shock propagating more rapidly into the low-density voids, instead of into the high-density filaments, of the cosmological density field 
(see also Fig. 7). While much of the swept-up mass is still outgoing at large velocities even after 50 Myr (pink and white), at this time 
metal-enriched gas is beginning to recoUapse into the host halo where it will likely form second-generation stars (see Fig. 6). 



heavy elements are mixed with the primordial gas, en- 
riching it to relatively high metallicity. 

Figure 5 shows the average metallicity to which the 
gas in the vicinity of the explosion site is enriched by 
the ejecta, after 1 Myr, 10 Myr and 70 Myr. As the 
blast wave propagates outward into the low-density IGM 
gas at greater and greater radii becomes enriched. As 
metals are carried out of the host halo, metallicities at the 
smallest radii fall. After 70 Myr, however, the average 
metallicity of the gas out to almost 10 kpc (physical) is 
enriched to of the order of 10~^ Zq, and the densest gas 
which is recoUapsing into the host halo is enriched to ~ 
0.05 Zq. As shown in Figure 6, by 70 Myr the gas at 
the center of the SN remnant, shown in Figure 7, has 
recollapsed to densities of n ~ 10^ cm~^, significantly 
more dense than the highest-density (n ~ 10 cm^^) gas 
that remained within the virial radius of the host halo 
~ 10 Myr after the SN. Not only is the gas recoUapsing 
after 70 Myr, but, as shown in the right panel of Fig. 6, 
a portion of the SN ejecta is entrained in this gas. If this 



ejecta is well-mixed with the dense primordial gas, then 
we expect its metallicity to be ~ 0.05 Zq, as indicated 
in Fig. 5. 

Gas enriched to such a high metallicity is predicted 
to readily fragment into low-mass stars (e.g., Bromm et 
al. 2001; Santoro & Shull 2006; Schneider et al. 2006), 
even in the presence of the elevated LW background ra- 
diation field expected in regions of the universe where 
SMSs form (e.g., Omukai et al. 2008). Therefore, we ex- 
pect that second-generation stars would form from this 
SMS SN-enriched gas, and that a large fraction of these 
would have masses low enough (< 0.8 Mq) that they 
could still be present in the Galaxy today. While the 
nucleosynthetic signature of very massive (i.e. 140-260 
Mq) Pop III SNe has yet to be uncovered in extremely 
metal-poor stars (Cayrel et al. 2004; Beers & Christlieb 
2005; Frebel et al 2005; Lai et al 2008; Joggerst et al 
2010; Joggerst & Whalen 2011), it may have been found 
in high-redshift damped Lyman alpha absorbers (Cooke 
et al. 2011), and a number of metal-poor stars in a re- 
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Fig. 5. — The spherically-averaged metaUicity Z of the gas (in 
units of the solar metaUicity Zq), as a function of distance r from 
the explosion site, at 1 Myr {blue dashed line), 10 Myr (green solid 
line) and 70 Myr (red dotted line). As the metal-enriched gas is 
blown out of the halo, the metaUicity in the interior region drops 
rapidly, while the IGM (outside the ~ 10^ pc virial radius of the 
host halo) is enriched with metals out to continually larger radii. 
After 70 Myr, the gas recoUapsing into the host halo (see also Fig. 
6) has a metaUicity of ~ 0.05 Zq. 
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Fig. 6. — The hydrogen nuclei number density n, as a function 
of the distance r from the explosion site, at 70 Myr after the SN. 
As in Fig. 2, the contours in the left panel denote the distribution 
of the gas, the mass fraction varying by an order of magnitude 
across contour lines. Individual SPH particles are shown in the 
right panel, with those constituting the metal-rich SN ejecta de- 
noted by orange triangles, as in Fig. 1. The highest-density gas is 
recoUapsing into the halo which hosted the SN, where it will likely 
form second-generation stars. While the majority of the metal-rich 
ejecta particles are in the low density intergalactic medium, the 
densest gas comprises some of them, suggesting that the second- 
generation stars will be highly metal-enriched. Many such stars 
would likely be long-lived and could still exhibit the chemical sig- 
nature of SMS SNe today. 



cent extension to the SEGUE survey have now been se- 
lected for spectroscopic fohowup on suspicion that they 
too may harbor this pattern (Ren et al. 2012). Fur- 
thermore, most stars forming in the ashes of very mas- 
sive primordial SNe were likely enriched to metallicities 
above those targeted by surveys of metal-poor stars to 
date (Karlsson et al. 2C)08), and our simulations predict 
that second-generation stars formed from gas enriched by 
SIVIS SNe will have metallicities above 10~^ Zq, which 
are also above this threshold. We conclude that, while 
SIVIS SNe are almost certainly rare events, their chemi- 
cal signature might be found in the ancient, metal-poor 
stars inhabiting our Galaxy today. 



4. DISCUSSION AND CONCLUSIONS 

We have carried out a multi-scale simulation of the 
explosion of a 55,000 Mq SMS, which, with an energy 
approaching 10^^ erg, is among the biggest explosions 
in the history of the universe. With our multi-code ap- 
proach, we have captured self-consistently, and for the 
first time, the essential radiation hydrodynamic features 
of the explosion at early times and the interaction of the 
blast wave with its host protogalaxy and with the IGM 
at later times. 

Although the atomic cooling halos expected to host the 
formation of SMSs are at least two orders of magnitude 
more massive than the DM halos in which the first pri- 
mordial stars are expected to form, we have found that 
these SMS SNe are energetic enough to completely evac- 
uate them of dense gas. The metal-enriched ejecta is 
dispersed well beyond the ~ 1 kpc virial radius of the 
host halo, out to ^ 5-10 kpc into the low-density IGM, 
after ^ 50 Myr. By this time, ^ 99 percent of the kinetic 
energy in the explosion has been radiated away; nonethe- 
less, the expansion of the shock into the IGM continues 
even at these late times, as the remaining kinetic en- 
ergy is still comparable to that expected for a ^ 200 Mq 
Pop III PSN. 

Because of the deep potential well of the halo and on- 
going accretion from filaments, after 70 Myr a fraction 
of the metal-enriched gas in the SN remnant, shown in 
Fig. 7, has recoUapsed to high densities. Given the rela- 
tively high metaUicity of this gas (~ 0.05 ZQ)0it is most 
likely that it will fragment vigorously and form a cluster 
of second-generation stars. Enriched to this metaUicity, 
any of these stars which survive to the present-day would 
exhibit metallicities much higher than the most metal- 
poor stars that have been found in surveys of the Galac- 
tic halo. Indeed, although SMS SN are likely rare events, 
it is possible that their chemical signature could be found 
in very old, relatively high-metallicity Pop II stars which 
inhabit the Milky Way today. As we expect SMS SNe 
to produce very little nickel-56 (Heger et al. 2013), their 
chemical signature may be distinguishable from those of 
most very massive (i.e. 140-260 Mq) Pop III explosions 
(PSNe), which can produce iron-group elements (Heger 
& Woosley 2002). Similar to these PSNe, however, SMS 
SNe would not make any s-process or r-process contribu- 
tions. 

Finding the chemical signature of SMS SNe in low- 
mass, long-lived stars (or in the IGM; see e.g., Cooke 
et al. 2011) would be one way of verifying that such 
exotic events occurred in the early universe. Another 
possibility is that these gargantuan explosions could be 
found in all-sky surveys such as those planned for the 
Wide-field Infrared Survey Telescope (WFIRST) and the 
Wide-field Imaging Surveyor for High-redshift (WISH), 
as shown recently by Whalen et al. (2012a)O Given 
that SMSs are expected to form in regions subjected to 
a large flux of LW radiation from nearby (within '^ 10 



^ Previous cosmological simulations of metal enrichment by mas- 
sive Pop III stars have shown second-generation star-forming gas 
to have typical metallicities of > 10~^ Zq (e.g., Wise & Abel 2008; 
Greif et al. 2010; Ritter et al. 2012; see also Wise et al. 2012). 

^ We note that neutrino emission from these explosions, pro- 
duced as discussed in e.g.. Yuan et al. (2011), might also be de- 
tectable. 
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Fig. 7. — The SMS SN remnant 70 Myr after the explosion (at z ~ 13), as seen on the sky from two orthogonal viewing angles. Due 
to the inhomogeneous cosmological density field into which the blast wave propagates, the remnant takes on an irregular shape which can 
subtend up to ~ 3 arcsec, making it many times larger than most galaxies at this early epoch. Indeed, it is likely that such remnants could 
envelop entire early star- forming galaxies (see Section 4). SMS SNe may be detected in all-sky surveys, such as those planned for WFIRST 
and WISH. At the center of this remnant, a cluster of metal-enriched second-generation stars is expected to form (see Fig. 6). 



kpc; Dijkstra et al. 2008; Agarwal et al. 2012) star- 
forming galaxies, detection of SMS SNe would pinpoint 
locations on the sky where rapidly-forming first galaxies 
could be found in follow-up observations by the James 
Webb Space Telescope. Indeed, these explosions are so 
large that the SN remnants they leave behind, as shown 
in Fig. 7, are likely to be many times larger than, and 
in fact are likely to envelop, any such neighboring star- 
forming galaxies. 

Could later stages of SMS SNe be detected by other 
means? They might appear in the radio at 21 cm. 
Meiksin & Whalen (2013) recently found that syn- 
chrotron emission from hypernovae in relatively dense 
media will be visible at 21 cm to existing radio facili- 
ties such as such as eVLA and eMERLIN in addition 
to the Square Kilometer Array (SKA). With their much 
higher energies and similar circumstellar densities, SMS 
SNe may be much brighter in the radio and detectable in 
all-sky surveys in spite of their small numbers. We are 
now calculating the radio signatures of SMS SNe in z '--^ 
15 protogalaxies. 

Likewise, Oh et al. (2003) and Whalen et al. (2008) 
examined the potential imprint of Pop III SNe on the 
cosmic microwave background (CMB) via the Sunyayev- 
Zeldovich (SZ) effect. They found that a population of 
140-260 Mq PSNe might impose excess power on the 
CMB on small scales, but with two caveats. First, the 
explosions must occur in large H ii regions because the 
SN shock must still be hot by the time it encloses a large 
volume of CMB photons. Explosions in dense media dis- 
sipate too much heat to have an appreciable SZ signature 
or upscatter many CMB photons. Second, although a 
population of Pop III SNe might collectively impose fea- 
tures on the CMB individual remnants achieve radii just 



below the current resolution of Atacama Cosmology Tele- 
scope or the South Pole Telescope. Our models explode 
in dense environments and at redshifts at which inverse 
Compton cooling losses are lower than for the first SNe, 
but they eventually reach radii that would allow them 
to be resolved by current instruments. We are currently 
evaluating the SZ signatures of SMS SNe. 
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